2/ci?-Friedel to Wigner oscillations in the one-dimensional Fermi gases under 

confinement 
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Density oscillations of confined one-dimensional Fermi gases of contact repulsive interactions in a 
continuous space are discussed within Bethe-ansatz based spin-density functional theory. The results 
are compared against the exact analytical and the exact diagonalization method. For unpolarized 
system, the number of the peaks in the density profiles is doubled signaling the crossover of the 
2fc_F-Friedel to 4fcF-Wigner oscillations (with fcp the Fermi wave vector). For both unpolarized and 
polarized systems, a threshold of the short-range interaction strength can be found where iV/-peak 
Wigner oscillations of 4fc_F wave vector appear in the density profile [Nf is the total particle number). 

PACS numbers: 03.75.Ss,71.15.Mb,71.10.Pm 



INTRODUCTION 



In higher dimensions a true crystal, named Wigner 
crystal, where the density is perfectly ordered, is possible 
when the repulsion between particles dominates over the 
kinetic energy in the case of small density. In one dimen- 
sion, quantum fluctuations prevent any type of true spin- 
symmetry breaking and long-range order in the ground 
state, and thus the crystal does not have perfect order 0- 
3]. However, in one dimension, the equilibrium positions 
of the electrons are assumed to be equally spaced and 
a solid like state can be produced when the correlation 
length exceeds the sample length [|| . The liquid-to-quasi- 
Wigncr crystal crossover when the density decreasing can 
be characterized by the onset of the 4/cp-peak in the 
static structure factor following from the Luttingcr liquid 
predictions 0, Q , with kp the Fermi wave vector (related 
to the fermion density n by kp = im/2), or by the Nf - 1 
high-contrast peaks in the pair correlation function Q, 
with Nf the particle number in the system. 

In the system where the constituent particles are not 
electrons but dilute atoms, the fundamental interaction is 
no longer of long-range Coulomb type but of short-range 
contact one, which is the case in the field of ultracold 
atomic systems. In this system, the ground-state den- 
sity profile displays at strong coupling well-pronounced 
Wigner oscillations with a l Akp periodicity'. However, 
the static structure factor or the pair correlation func- 
tion does not show any signature of Wigner-molecule- 
type correlations like that in the system of long-range 
interaction Q, which means Wigner- type molecules are 
absent in the system of short-range interaction Q . In this 
paper, we focus only on the :: 2kp -*■ Akp crossover" in the 
density profiles: for vanishing and weak interactions the 
ground state is liquid-like with 2fcF-Friedel oscillations, 
whereas at strong repulsions a periodicity of 4fcp-Wigner 
oscillation emerges. 



For the one-dimensional (ID) harmonically trapped 
system of spin-balanced or imbalanced fermionic atomic 
gases, many numerical schemes have been used in study- 
ing lots of interesting quantum effects pj-[l4l] . A density- 
functional scheme using as reference the Bethe-ansatz 
solvable homogeneous system was proposed to calculate 
the ground-state properties [HI, [Tg. However, this ap- 
proach was found to fail at the strong-coupling limit [161 ] . 
More precisely, the exchange-correlation potential pro- 
posed there is not able to describe the 2kp -»■ Akp 
crossover, which, as noted above, is expected to occur 
upon increasing the strength of the repulsive interactions 
between antiparallel-spin particles. 

A simple functional is proposed in Ref. 
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particle system that embodies this crossover and is capa- 
ble of describing inhomogeneous Luttinger liquids with 
strong repulsions. The main idea is to capture the ten- 
dency to Akp- Wigner oscillations by adding an infinitesi- 
mal spin-symmetry-breaking magnetic field to the Hamil- 
tonian and by resorting to spin-density functional theory 
(SDFT) based on the spin-polarized Bethe-ansatz equa- 
tions mm- 

For fermions loaded in a ID optical lattice, the system 
can be modeled by a ID Fermi-Hubbard model. Friedel 
oscillations of wavelength 2kp and Akp component are 
analyzed in this model with open boundary condition by 
both the Bethe-ansatz based SDFT, the exact diagonal- 
ization [20j . bosonization, and numerical density matrix 
renormalization group calculations [2lTj23| . It was found 
that the 2kp -* Akp crossover takes place at U = At. For 
an unpolarized system, if a tiny magnetization which is 
proportional to the difference of the spin-up and spin- 
down densities is induced by hand into the exchange- 
correlation potential, the Akp -oscillation appears in the 
density profile at the expense of violating the Lieb-Mattis 
theorem, while for a polarized system, SDFT can repro- 
duce both the 2kp and Akp oscillations without violating 
the Lieb-Mattis theorem [2(|. It is found that the Akp 
density oscillations is strongly related to the large Fermi 
sea induced by the strong repulsion between different spin 
species pll. [23], 

Some ID systems, for example, the ID Gaudin-Yang 
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gas in a homogeneous continuum space, or the ID Lieb- 
Wu model of particles in a discrete lattice, i.e., ID Fermi- 
Hubbard model, where the fcrmions interact through at- 
tractive or repulsive short-range interactions, are exactly 
solvable. However, in the inhomogeneous case only some 
special systems can be solved exactly. For example, the 
ID two-component Fermi gas in an infinite potential well 
can be solved exactly by the Bethe-ansatz techni ques , 
where the ground state properties are investigated [24| j . 
The density profiles are evaluated from weak to very 
strong interactions. In the strongly interacting limit, 
the density profile tends to that of the spinless fermions 
though the momentum distributions differ. An exact an- 
alytical solution of ID spin-1/2 fermions with infinite 
repulsion for arbitrary confining potential is obtained 
by the combination of Girardeau's hard-core contacting 
boundary condition and group theoretical method [25| . 
For the spin-imbalanced case they found that the spin- 
up and spin-down density distributions show alternative 
peaks avoiding overlapping together and the parity sym- 
metry. 

In this paper, we will employ the Bethe-ansatz based 
SDFT to investigate the ID inhomogeneous system with 
general interaction strengths. Our results are compared 
against the above mentioned exact analytic and the exact 
numerical diagonalization method. The Wigner oscilla- 
tions while increasing coupling strength for both unpo- 
larized and polarized cases will be studied. 

The rest of the paper is organized as follows. In Sec. 
II, we present the model and the method. The numerical 
results are followed in Sec. III. In Sec. IV, we conclude 
what we found. 
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FIG. 1: (Color online) The ground-state density distributions 
as a function of z/au for a balanced Fermi gas with AT-f = 
Ni = 5 confined in a harmonic trap for various interaction 
strengths. The crossover from 2k p- to 4fci?-oscillations with 
increasing interaction strength A is showed by the doubling 
of the crests in the density profile. The dominance of the 
4fci?-oscillations over 2fci?-ones can be better evidenced by the 
density derivative dn(z)/dz. In the inset, the derivative of the 
density with respect to z is shown for A = 1,6, and 8 as an 
indicator to examine the 2&f -* 4&f crossover. 



Hamiltonian ([I]) can be shown to be governed by the di- 
mensionless coupling parameter 



II. MODEL AND THE METHOD 



A 



Q,|| 



(2) 



In this work we consider a ID system of two-component 
Fermi gas with a contact repulsive interaction of strength 
5id [26j j . The two fermion species are taken to have the 
same mass m, with Nf spin-up species and N± spin-down 
ones. The total number of particles is Nf = iVf + iVj,. The 
system is described by the Hamiltonian, 



2m U dz\ 



H = ~— E — + SID E E S(Zi ~Zj) + Kxt , (1) 

i=lj=l 



where contact interactions between parallel-spin particles 
are suppressed by the Pauli exclusion principle and V ex t 
the confining potential. We consider that the atoms are 
either confined in a ID box of length L (V^ xt = inside the 
box and V cx t = +°° elsewhere) or subjected to a strongly 
anisotropic harmonic potential V4 x t = Zj=i 14xt(^i) = 
Y,i=i z fi characterized by angular frequencies us ± 
and wii in the radial and axial directions with wi « lo l . 



For the system of harmonic trap, we choose the harmonic- 
oscillator length a || = yjh/(mujj) as unit of length and the 
harmonic-oscillator quantum hu>\\ as unit of energy. The 



For Vcxt = the Hamiltonian (fTJ) reduces to the ho- 
mogeneous Gaudin-Yang model solvable by means of 
the Bethe-ansatz technique [27j and determined by the 
linear density n = N/L, by the spin polarization £ = 
(iV-f - N±)/N, and by the interaction strength g^. 

The ground-state spin-density distributions n a (z) can 
be calculated within SDFT by solving sclf-consistently 
the Kohn-Sham (KS) equation, 



r<Pa,<r(z) (3) 



together with the closure 



N a 

E \<p<xAz 



Here, V^[n a ]{z) 
is the spin-dependent effective KS potential 



]{z) + V^\n a \{z) 



(4) 

+ V cxt (z) 
The first 

term in V^g is the mean-field term = giT>n s (z), 

while the second term is the exchange-correlation poten- 
tial defined as the functional derivative of the exchange- 
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FIG. 2: (Color online) The ground-state density distribu- 
tions (in arbitrary length -1 unit) as a function of z (in ar- 
bitrary length unit) for a polarized Fermi gas with Nf = 2 and 
ATj = 1 confined in a box of length L = 1 for various inter- 
action strengths. The SDFT results are plotted in symbols 
and the exact Bethe-ansatz ones in (colored) lines. The thin 
(black) solid line represents the result in the infinite repulsive 
interaction limit. 



correlation energy E xc [na-] evaluated at the ground-state 
density profile, V^c = 6E xc [n a ]/6n,j(z)\ GS . 

To obtain the density distribution n a {z) from Eqs. © 
and ©, approximations to the exchange-correlation 
functions E KC are unavoidable. The local-spin-density 
approximation (LSDA) is known to provide an excellent 
description of the ground-state properties of a variety of 
inhomogeneous systems (28j . In the following we em- 
ploy a Bethe-ansatz-based LSDA (BALSDA) functional 
for the exchange-correlation potential, 

E xc [n a ] -» £ X L C SDA K] (5) 
= J dzn(z) £xc° m KC,3iD)L n(z) ^ c(z) , 

where the exchange-correlation energy e^° m of the homo- 
geneous Gaudin-Yang model is defined by 

e x? m ( n >C,5iD) = £gs(™,C.9id) - n(n,() 

- £h(™,C,3id)- (6) 

Here ecsin, C;<7id) is the ground-state energy of the ex- 
act Bethe-ansatz solution of the model, en(n, C,<?id) = 
jgiDn 2 (l - C 2 ) is the mean-field energy, and ft(rt,0 = 
7r 2 fc 2 n 2 (l + 3C 2 )/24m is the noninteracting kinetic en- 
ergy [27|. 



III. NUMERICAL RESULTS 

In this section, we present our numerical results by 
employing the Bethe-ansatz based SDFT to investigate 
the ID inhomogeneous system with broad interaction 
strengths. Our results are compared against the exact 
numerical diagonalization method while this method is 
available when the system is small. For the polarized 
system in a box, wc compare our results with the ana- 
lytical ones which are available in this case by the exact 
Bethe-ansatz method. 



A. Unpolarized system confined in a trap 

In Fig. [TJ we study an unpolarized Fermi gas with the 
number of total particles Nf = 10 confined in a harmonic 
trap. The ground-state density distributions is calculated 
for various interaction strengths within SDFT. For the 
unpolarized system, the KS equations have to be solved 
sclf-consistcntly by taking one of the following methods 
to obtain the correct 2kp -»■ 4/cf crossover, (i) One can 
cither add an infinitesimal spin-symmetry-brcaking mag- 
netic field. The field strength is infinitesimal small and 
reduced to zero during the self-consistent solution of the 
KS equations, (ii) Alternatively one can solve the KS 
equations from an initial guess which has slightly broken 
spin symmetry [i.e. an initial guess with n^(z) + 
without using any B field, (iii) The last possibility is to 
break the symmetries of the density distributions just one 
time by shifting the center of the densities during the 
self-consistent cycle [for example, let n-f(z) = n^(z + b) 
with b a nonzero value]. When A is small, the above 
mentioned methods restore the results from the Bethe- 
ansatz based local density approximation (B ALDA) pj|] . 
However when A is large, the final result for the total 
density n(z) from the BALSDA is very different from 
that resulting from the B ALDA theory [16[ . We want to 
mention here that, without using the above mentioned 
symmetry-breaking procedure, the BALSDA theory pro- 
duces the same results as that of BALDA. That is, for 
the BALSDA theory, to achieve a correct shape for the 
density profile, one has to use one of the (i)-(iii) proce- 
dure This method has long been used in the liter- 
ature in discussing the ground-state of the H2 molecule 
at any bond length. The symmetry-breaking choice typi- 
cally produce a spin-unpolarizcd density out to a critical 
bond length by localizing an up-spin electron on one nu- 
cleus and a down-spin electron on the other. While this 
leads to an incorrect spin-density, but it leads to the cor- 
rect dissociation limit for the energy and the correct total 
density [H. 

At the weak interaction regime, the density is charac- 
terized by 2/c^-Friedel oscillations with Nf/2 peaks. As 
the interaction strength increases to the strongly inter- 
acting limit, the number of peaks in the density profile is 
doubled. The atoms are distributed in order to minimize 
the interparticle repulsion energy. This leads to the emer- 
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FIG. 3: (Color online) The ground-state density distributions as a function of z/au for a polarized Fermi gas with N\ = 2 and 
JVj = 1 confined in a harmonic trap ranges from weak to strong interaction strengths: (a) A = 1; (b) A = 5; (c) A = 8; (d) A = 15. 
The SDFT results are plotted in symbols (°, », and A, respectively, for n, rif, and nj ) and the exact diagonalization ones in 
lines. 
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FIG. 4: (Color online) The 3D plot of the ground-state density 
distributions as a function of the dimensionless interaction 
strength A and z/a^ for a polarized Fermi gas with N-\ = 2 and 
Ni = 1 confined in a harmonic trap. The results are obtained 
based on the SDFT. 



gence of Nf crests in the density profile characterized by 
strong correlations and the 4/cF-Wigner oscillations. In 
this case the atoms behave like a fully polarized (or non- 



interacting spinless) fermions as in the Tonks-Girardeau 
gas Numerically we find that a crossover from the 
2fcp-Friedel oscillations of Nf/2 peaks to the 4fcp--Wigncr 
oscillations happens around A = 6.2 accompanied of Nf 
peaks. The crossover is smooth, see, for example, the 
extensive density- matrix renormalization group study by 
Soffing et al. [2l|, [22|]. In the inset, a derivative of the 
density with respect to z/aj is used as an indicator to 
signal the 2kp -> 4fcp crossover. We find more clearly 
that Nf peaks appear after the threshold value. 

In order to evaluate the credibility of the present 
method, we compare our results in the following against 
exact methods for a polarized system in a box (in IIII Bp 
and a harmonically confined system of small number of 
Fermi atoms (in IIII Cjl . 



B. Polarized system in a box 

In this section we consider a polarized Fermi gas with 
Nf up-spin fermions and N± down-spin fermions con- 
fined in a box of length L, which is exactly solvable by 
the Bethe-ansatz method [241] . In Fig. [2J we illustrate 
the ground-state density distributions as a function of 



FIG. 5: (Color online) The 3D plot of the ground-state density 
distributions as a function of the dimensionless interaction 
strength A and z/a« for a polarized Fermi gas with Ni = 3 and 
Ni = 2 confined in a harmonic trap. The results are obtained 
based on the SDFT. 



FIG. 7: (Color online) The 2D contour plot of the ground- 
state density distributions as a function of the dimensionless 
interaction strength A and z/a^ for a polarized Fermi gas with 
ATf = 34 and N± = 6 confined in a harmonic trap. The results 
are obtained based on the SDFT. 




FIG. 6: (Color online) The 2D contour plot of the ground- 
state density distributions as a function of the dimensionless 
interaction strength A and z\a\\ for a polarized Fermi gas with 
N\ = 22 and JVi = 18 confined in a harmonic trap. The results 
are obtained based on the SDFT. 



z based on the Bethe-ansatz based SDFT for N\ = 2 
and JVj. = 1. L is chosen as units of length. The in- 
teraction strength sets the relative importance of the ki- 
netic energy and the interaction. For weak interaction 
strength, 2 peaks appear in the density. When inter- 
action strength becomes stronger (around gm = 15), 3 
peaks emerge. At even stronger interaction strength, ex- 
change is suppressed and the atoms act as though they 
were nonintcracting and spinless in analogy to the case 
studied for bosons with an infinite contact repulsion [3Cjj . 



For gio = 1000, the infinite repulsive interaction limit is 
reached, which is compared to the spinless fermions re- 
sults. In general, when the interaction strength is large 
enough, the density n(z) distributes between Nf max- 
ima, which is different from the confined ID electrons 
of long-range Coulomb interactions where the density 
vanishes almost completely between well separated max- 
ima [U [HI]. We find, in almost all the range of the 
interaction strength, SDFT gives good results comparing 
against the exact Bethe-ansatz ones. The good agree- 
ment may be that the exact exchange-correlation energy 
of the homogeneous system has almost the same form as 
that of the inhomogeneous (hard- wall) one which is also 
integrable. 

C. Polarized system confined in a trap 

Here we report our studies on the ground-state den- 
sity distributions for a polarized Fermi gas confined in a 
harmonic trap. The SDFT results are compared against 
the ones by the exact diagonalization method [25| . In 
Fig. [3] we plot the density profiles for the system with 
JV t = 2 and JVj, = 1 for A = 1,5,8, and 15. We find that 
for A < 2, both the total ground-state densities n{z) and 
the spin-resolved densities show very good agreement be- 
tween the SDFT results and the exact diagonalization 
ones. When A > 2, SDFT still keeps the overall shape 
of the total ground-state densities well. However, the 
respective amplitudes of the spin-up and spin-down den- 
sity are exaggerated by the SDFT. Even for this weak- 
ness, the SDFT can help us qualitatively understand the 
2kp -> 4fcp crossover. For weak interaction strength, the 
interaction term gives smaller contribution comparing to 
the kinetic term, such that n(z) is dominated by the 
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lowest occupied single-particle states. For the system 
we studied here, one flat peak is found (See Fig. E3a) 
and (b)). For strong interaction strength A > 3.8 [33|], a 
structure consisting of Nf = 3 peaks emerges (See Fig. 
[3^c) and (d)), that is, n(z) is signaled in finite region by 
Nf = 3 maxima. For even stronger interaction strength 
A > 15, the total atom distribution can be described by 
spinlcss fermions and the spins tend to decouple with 
exponentially small exchange interactions. The system 
degenerates with ferromagnetic and anti-ferromagnetic 
states [25j . 

In Figs. [Hand 03 the 3D plots of the ground-state den- 
sity distributions are shown as a function of the dimen- 
sionless interaction strength A and z/au. The transition 
from 2fcp-Friedel oscillations to 4fcp-Wigner oscillations 
can be judged by the change in the peaks in the total 
density distributions, which vary from 1 to 3 in Fig. [4] 
happened at A > 3.8, and from 3 to 5 in Fig. [5] happened 
at A> 4.7. 

At strongly repulsive limit, the polarized system shows 
the same properties as those of the unpolarized one, 
where the atoms become impenetrable and behave as 
they were noninteracting and spinlcss. 

DFT provides a method through which we can access 
larger systems beyond the reach of QMC and DMRG 
simulations. In the following we study the polarized sys- 
tem confined in a trap with particle number of Nf = 40, 
which is not accessible by means of the present exact di- 
agonalization method. Contour plots arc shown in Figs. 
[6] and [3 respectively, for p = 0.1 and 0.7. A clear signa- 
ture of 2kp -+ 4kp crossover is signaled by the emergence 
of Nf peaks at the strong interaction strength. 

From the above discussion, we find that, in spite of 
the short-range nature of the interaction, density profile 
shows pronounced well-separated Nf peaks at strong re- 
pulsive interactions. However, the densities between the 
peaks do not drop to zero, which is different from the ID 
system of long-range Coulomb interaction, where they 
drop almost to zero in some sense resembling a charge 
density wave and a quasi- Wigner molecule. 

IV. CONCLUSIONS 

In this work, we study a ID system of two-component 
Fermi gases with a contact repulsive interaction. The 



unpolarized or polarized atoms are either confined in a 
ID box or subjected to a strongly anisotropic harmonic 
potential. The total and spin-resolved densities are cal- 
culated within the Bethe-ansatz based spin-density func- 
tional theory. The results are compared to those from the 
exact Bethe-ansatz method and the exact diagonalization 
scheme. We found that the Bethe-ansatz based SDFT 
gives a faithful description for the ground state prop- 
erties at weak interaction strength but exaggerates the 
amplitude of the density distributions especially for the 
corresponding spin-up and spin-down densities at strong 
interaction strength. In spite of the short-range nature 
of the interaction, a pronounced iV/-peak structure is 
observed in the density profile in the strongly interacting 
regime. 



For unpolarized systems, the density is characterized 
by 2/cF-Friedel oscillations with Nf/2 peaks at the weak 
interaction regime. As the interaction strength increases, 
the number of peaks in the density profile is doubled. 
This leads to the emergence of 4k f- Wigner oscillations of 
Nf peaks in the density profile. Numerically we found a 
threshold interaction strength A c at which the 2kp ->• 4kp 
transformation happens. For polarized systems, at weak 
interaction strength, the peaks in the density profile de- 
pend on the system parameters. However, when increas- 
ing the interaction strength, a threshold value can be 
found where A^/-peak oscillations appear in the density 
profile signaling the 4fc^-Wigner oscillations. At strongly 
repulsive limit, the total energy and the total density pro- 
file of the system own the same properties, as those of 
the fully unpolarized one. There, the atoms become im- 
penetrable and behave as they were noninteracting and 
spinlcss. 
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